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When an epidemic spreads in a population, individuals may adaptively change the structure of 
their social contact network to reduce risk of infection. Here we study the spread of an epidemic 
on an adaptive network with community structure. We model the effect of two communities with 
different average degrees. The disease model is susceptible- infected-susceptible (SIS), and adaptation 
is rewiring of links between susceptibles and infectives. The bifurcation structure is obtained, and 
a mean field model is developed that accurately predicts the steady state behavior of the system. 
We show that an epidemic can alter the community structure. 
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I. INTRODUCTION 

In recent years, networks have been widely used in mod- 
eling a variety of social, technological, and biological sys- 
tems [H-Q . One major application is modeling the spread 
of an epidemic on a social network (e.g., [4-8]). In these 
models, typically the network structure is assumed static, 
while the infection status of the nodes changes dynami- 
cally. 

During an epidemic, people may tend to avoid social con- 
nections with infected individuals d, This system 
can be considered an adaptive network, where the node 
dynamics affects the network topology, which then affects 
future changes in node status [ll|- Epidemic spreading 
in adaptive network models with avoidance behavior has 
been studied previously || IIlT - fl5j . with avoidance fre- 
quently implemented via susceptible nodes rewiring their 
links away from infected neighbors and towards other 
non-infected nodes. Changes in bifurcation structure 
have been observed, including the existence of bistable 
regimes with endemic and disease free states both stable. 
One of the most important features of a social network is 
community structure [l6l |. The strength of the commu- 
nity structure can be quantified by using a modularity 
measure [l6| . which for a random network will be close 
to zero, and will be close to one for a strong community 
structure. 

Studies of epidemics with community structure have fo- 
cused mainly on static network geometries, including 
scale- free p!7H19| . small- world [20| and random networks 
[2ll ] . It has been found that community structure can ei- 
ther decrease [l8| or increase [2l[ infection prevalence, de- 
pending on details of the model. Further, epidemics can 
synchronize across communities if there are sufficient con- 
nections between communities [l7l . |20| . In a dynamic but 
not adaptive example, communities of mobile agents were 
studied [13, HI] , and dynamic hopping of agents between 
communities was able to produce sustained infection in 
communities that were below the epidemic threshold if 
other communities were above threshold. 
Adaptive networks with community structure have been 
studied only rarely. In [24j . the authors considered an 
adaptive scale-free network with community structure in 
which neighbors of an infected node can move to other 
communities with a certain probability. Infection levels 



were reduced compared to the case without adaptation, 
but the adaptation mechanism did not preserve the com- 
munity structure as measured by modularity. In [251 ] , the 
authors introduced a model very similar in structure to 
the one we will consider here. However, their focus was 
to study an adaptive epidemic system with two types of 
agents. They varied the within- type and cross- type link 
rewiring rates and infection rates and determined their 
effects on the size of the bistability region. For certain 
parameter choices the endemic steady state would have 
a community structure, but the resulting structure was 
not characterized in this study. 

In this paper, we extend the adaptive susceptiblc- 
infected-susccptiblc (SIS) model of [9( to a network with 
two communities. In contrast to previous studies [23.[25j. 
we allow the communities to have different average de- 
gree. We define rewiring rules such that the community 
structure is preserved if links between susceptibles and in- 
fectives are uniformly distributed. We directly simulate 
the stochastic network system and derive a lower dimen- 
sional mean field, based on a moment closure approxima- 
tion, that accurately predicts the bifurcation structure of 
the full system. In Section [UJ we define the model and 
introduce the mean field equations. Results in the ab- 
sence of adaptation are presented in Subsection IIII Al In 
Subsection IIII Fit we show the effects of adaptation on 
the bifurcation structure and on the network geometry. 
Section [TV] concludes. 



II. MODEL 

We study a susceptible-infective-susceptible (SIS) model 
on an adaptive network having two communities. The 
communities are labeled A and B and consist of N a and 
Nf, nodes, respectively. Here N a = Nf, — 5000. We 
use two probability parameters to generate an initial net- 
work with two communities by creating links. Parameter 
d € [0, 1] determines the asymmetry in average degree of 
the communities, and / € [0, 1] determines the number 
of links between communities. Links are created as fol- 
lows. With probability d we choose a node among the 
N a nodes in community A (otherwise choosing a node 
in community B), and with probability / its neighbor is 
chosen at random from the opposite community as the 
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first node (otherwise choosing from the same commu- 
nity). This process is repeated until a total of K links 
are created. Self links and multiple links arc disallowed. 
Then a fraction ( 1 — / ) d of the links are A A , ( 1 — /) (1 — d) 
are BB, and / are AB. The average degrees in commu- 
nities A and B are (k a ) = [2(1 - f)d + f] K/N a and 
(k b ) = [2(1- f)(l-d)+ f]K/N b , respectively. Thus 
the communities are symmetric when d = 0.5. We will 
focus here on the case d > 0.5, so community A will have 
higher connectivity than B. The total number of cross 
links between communities is fK. 

We define node dynamics as in 0. A susceptible (S) 
node becomes infected with rate pNi n f, where Ni„f is 
the number of infected neighbors the node has and p is 
the infection rate. An infected (I) node recovers with 
recovery rate r. One of these rates can be eliminated by 
rcscaling time, so it is sufficient to treat r as fixed. We 
fix r = 0.002 throughout the paper as in previous studies 

mm ■ 

Network adaptation in the form of avoidance behavior 
is introduced by allowing susceptible-infected links to 
rewire with rate w to susceptible-susceptible links, as 
in Q. However, the rewiring must be adjusted to re- 
tain the desired community structure. This is done by 
choosing the susceptible node's new neighbor from one 
or the other community with appropriate probabilities. 
An S node having an infected neighbor rewires to an S 
node in the same community as itself with probability 
a if the S node is in community A and with probability 
(3 if the S node is in community B. Otherwise a neigh- 
bor in the other community is selected. In order to re- 
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tain the community structure, we set a = 2 (i-f)d+f 

P = 2 {i-~fj(i-~d)+f ■ This choice is made so that if ran- 
domly selected links rewire, then the flux from AA links 
to AB links, (1 — f)d(l — a), equals the flux from AB 
links to A A links, ^a, and likewise for balance of fluxes 
between BB and AB links. Therefore, if SI links occur at 
random anywhere in the network, this rewiring strategy 
will on average keep the community structure specified 
above by d and /. 

We simulate our model using Gillespie's method [26| for 
N = 10 4 nodes and K = 10 5 links The initial condi- 
tion is either the final state of a previous run or a random 
two-community network constructed as described above 
in which a fraction of the nodes have been randomly in- 
fected. 

As in 0, [l3| , we derive mean field equations for the evo- 
lution of the nodes and links. Px denotes the probability 
of nodes to be in state X, where X is susceptible in com- 
munity A or B (S a or S b ) or infected in A or B (I a or It). 
Pxy denotes the probability that a randomly selected 
link connects a node in state A to a node in state Y. We 
obtain the following equations for the node dynamics: 



pK 

-rPl a + — (Psala + Psalb) 
^ * a 

pK 



(1) 



Pit, = -rP Ib + ?—(Psbia + Psbib) (2) 



Because nodes are neither created nor destroyed and do 
not change their community assignment, the equations 
for susceptiblcs in community A and B can be found from 
node conservation. 

The evolution of the links depends on three point terms. 
As in 0, [TH, [TU we use a moment closure assumption 
to close the system, assuming Pxyz ~ PxY p ^ YZ , where 
Pxyz is the fraction of three point terms. After applying 
the moment closure, the link equations are 
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These mean field e quat ions are a special case of the gen- 
eral mean field in [25[ for appropriate choices of their 
rewiring and infection parameters. 

Since the total number of links is fixed, we have an 11 
dimensional system in the adaptive network case by elim- 
inating one of the link equations. On the other hand, we 
have a 9 dimensional system in the static network case 
because the numbers of AA, AB, and BB links are each 
fixed. These equations can be integrated with standard 
numerical integration methods. Also, we tracked their 
steady states using a continuation package [27j . 



III. RESULTS 
A. Static network 

We first consider the effect of having two communities 
with different average connectivities in a static network 
(w = 0). We obtained the bifurcation structure (Fig- 
ure [JJ) as follows. We used the XPPAUT free software 
package [27j to locate the stable and unstable equilib- 
rium solutions of mean field equations. To obtain the 
steady states of the full system, we generated an initial 
random network with community structure in which 50% 
of the nodes were infected. To locate the upper branch 
(endemic state), the system was run to steady state for 
a high infection rate p, and p was decreased gradually 
using the final state of each run as an initial state for the 
next run. For each p, we run the system up to 5 x 10 4 
time units and then averaged the steady state over 500 
samples where there are 10 3 events between each sample. 
To locate the lower branch (disease- free state), we gen- 
erated a random network with community structure in 
which 2.5% of the nodes were infected. The system was 
simulated for 5 x 10 4 time units, and five runs were done 
for each p value. If the infected fraction went to zero in 
any of the five runs, the disease-free state was considered 
stable. As shown in Figure [TJ the mean field equations 
and the full system are in good agreement. 
In a static network without community structure, the 
disease free state (DFS) loses stability at a critical infec- 
tion rate p* where the system undergoes a transcritical 
bifurcation. This threshold infection rate depends on the 
average degree of the network [9j . Figure QJ, superim- 
poses the bifurcation diagrams of two single-community 
networks with different average degrees. The epidemic 
threshold for community A (p*) and community B (p^) 
are significantly different. When the two networks are 
loosely connected (/ = 10~ 4 , Figure [T]d and blowup in 
Figure [T]i) , the combined system has a single threshold 
infection rate, which is approximately p* and much lower 
than pi in the disconnected case. When the infection 




FIG. 1: Bifurcation diagram for infected fraction in a static 
network (w — 0) as a function of infection rate p for different 
fractions of cross links /. Black: community A; gray: commu- 
nity B. Solid curves: mean field solutions (stable branches); 
dashed curves: mean field solutions (unstable branches); cir- 
cles: Monte Carlo simulations. Average degrees (k a ) w 28 and 
(k b ) « 12. The plots correspond to N a = 5xl0 3 , N b = 5xl0 3 , 
K = 10 6 , d = 0.7 and r = 0.002. In (d), in order to obtain 
better statistics, we averaged over 100 network realizations, 
each of which is similar to (b). 



rate is between p* and pi , the fraction of infecteds in B 
is close to zero (Figure QJI) and stochastic reintroduction 
of infection from A to B is observed. However, when the 
two communities are strongly connected (/ = 10~ 4 , Fig- 
ure [TJ;), they behave similarly in that both communities 
have significant infection levels for the same parameter 
values. 

Although a system of two connected communities has an 
infection threshold at a single bifurcation point, we wish 
to distinguish between the cases in Figure [1}d,c, where 
significant infection spread occurs in the low degree com- 
munity (B) if it has sufficient links to the high degree 
community (A), while the infection in community B is 
very small if the number of cross links is low. To quan- 
tify this, we define effective threshold infection rates p e a: 
p\ for each community as follows. While sweeping the in- 
fection rate p from higher to lower values, the first p value 
at which the fraction of infecteds at the steady state is 
lower than e = 10~ 3 is considered as the effective thresh- 
old infection rate for that community. Figure [2] shows the 
effective threshold infection rates versus cross link frac- 
tion /. When / < 10~ 2 , the effective thresholds in the 
two communities become noticeably different. 
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FIG. 2: Effective threshold infection rate from mean field vs 
fraction of cross links /. Black curve: community A; gray 
curve: community B. The infection rate p at which the frac- 
tion of infecteds is less than 10~ 3 is considered as the effective 
threshold rate. N a = 5 x 10 3 , N b = 5 x 10 3 , K = 10 6 , d = 0.7 
and r = 0.002. 



B. Adaptive network 

We now move to systems with nonzero rewiring rates. 
The bifurcation structure for the adaptive network case 
was determined as for static networks except for the fol- 
lowing modification. In our model, network adaptation 
does not occur in the absence of infection. This means 
that the DFS of Equations (|TT fT2|) is not isolated, be- 
cause any disease free combination of AA, AB, and BB 
links is a steady state. Because of the non-isolated fixed 
points, the stability of the disease free branch could not 
be determined using continuation packages. Instead, we 
calculated numerically the eigenvalues of the Jacobian 
evaluated at the DFS for the initial network geometry 
which is described by / and d. 

In the absence of community structure (/ = 0), the DFS 
loses stability at a critical infection rate p* , where the un- 
stable endemic branch and the stable disease free branch 
intersect at a transcritical bifurcation (Figure The 
epidemic threshold p* is inversely proportional to the av- 
erage degree of the network and can be found analytically 
from the Jacobian of the mean-field equations for a single 
network (see Appendix). 

For a network having two loosely connected heterogenous 
communities (/ = 10~ 4 , d = 0.7), p* is very close to that 
of a single network having the same average degree as 
community A (Figure [3^, b). This is expected because 
when the infection in community A starts to spread, a 
small value of / will not be enough to stop infection 
spreading in community A. In contrast, for / = 10" 1 , p* 
is larger than that of a single network having the same 
average degree as community A. In this case, the con- 
nection is stronger and the infection starting to spread 
in community A can be suppressed by the connection 
to a community where no infection is observed. As we 
increase /, the critical value of p approaches that of a 
single network with average degree (k) = 20, which is 
the average degree in the entire system. 
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FIG. 3: Bifurcation diagram for infected fraction in an adap- 
tive network as a function of infection rate p for different 
fractions of cross links /. Black: community A; gray: commu- 
nity B. Solid curves: mean field solutions (stable branches); 
dashed curves: mean field solutions (unstable branches); cir- 
cles: Monte Carlo simulations. Average degrees (k a ) ~ 28 and 
(k b ) w 12. The plots correspond to N a = 5x 10 3 , N b = 5x 10 3 , 
K = 10 6 , d = 0.7, w = 0.04 and r = 0.002. The simulations 
are done similarly as in Figure [U but here we run the sys- 
tem for 5 x 10 6 time units in order to approach the endemic 
branch. 



In an adaptive network without community structure, 
bistability can occur for a range of rewiring rates Q . We 
focus on the rewiring rate w — 0.04, for which the en- 
demic branch loses stability at a critical infection rate 
where the system undergoes a saddle-node (SN) bifur- 
cation. The location of this SN bifurcation depends on 
the average degree of the network. Figure [3^ superim- 
poses the bifurcation diagrams of two single-community 
networks with different average degrees. 
In our model, for a network having two strongly con- 
nected communities (/ = 10 _1 ), the endemic state loses 
stability at a SN bifurcation point. However, for the 
loosely connected case (/ = 10~ 4 ), the endemic branch 
loses stability at a Hopf bifurcation (HB) point. As 
with the epidemic threshold (transcritical bifurcation) in 
static networks, the location of the SN bifurcation and 
the HB in the adaptive network is governed primarily by 
the high degree community. However, in contrast with 
the static case, small cross link fraction / is not asso- 



5 




time 



x 10 



FIG. 4: Time series from Monte Carlo simulation of full sys- 
tem, (a) Infected node fractions. Black curve: fraction of 
infecteds in community A; gray curve: fraction of infecteds 
in community B. (b) Fraction of links. Black curve: frac- 
tion of AA type links; gray curve: fraction of BB type links. 
Initial condition is as described in the text. N a = 5 x 10 3 , 
N b = 5 x 10 3 , K = 10 6 , / = 1(T 4 , d = 0.7, w = 0.04, 
r = 0.002 and p = 0.0014. 



ciated with low steady state infection levels in the low 
degree community. As will be seen later in this section, 
the state with high infection in one community and low 
in the other community is not a steady state due to the 
network adaptation. Instead, infection levels are simi- 
lar in both communities even if weakly connected. Both 
communities continue to exhibit high infection levels as 
the number of cross links is further increased (FigureJSJ;). 
By looking at the time series of the mean field equations, 
we observed a stable periodic solution for a very small 
range of p values near the HB point. However, we do not 
see any periodic behavior in the full system due to the 
narrow range of p values. 

As the heterogeneity in the network increases by increas- 
ing d, the Hopf bifurcation point also increases. Further- 
more, the epidemic threshold decreases because of a much 
higher average degree in community A. This narrows the 
region where both disease free and endemic branches are 
stable. In particular, for d = 0.9, there is no p value 
where both the endemic steady state and DFS are sta- 
ble. For d = 0.9 We observed periodic solutions with a 
very long period for p values smaller than the HB point 
in both mean field and the full system. 
To motivate the absence of the state seen in static net- 
works with high infection in one community and low in- 
fection in the other, we consider a long time series start- 
ing initially with a network where / = 10~ 4 and d = 0.7 
generated as described above. As seen in Figure HI the 
fraction of AA and BB type links change with time. The 



adaptation rules have been chosen so that if SI links are 
distributed uniformly throughout the network, the com- 
munity structure will be preserved. However, with high 
initial infection levels in A and low in B, there are more 
SI links among the AA and AB links and fewer among the 
BB links. This leads to a net flux of link types from AA to 
AB to BB. (The fraction of AB links (not shown) remains 
relatively constant at low levels throughout.) Eventually 
the average degree in the B community exceeds that in 
the A community and there is an incursion of infection 
from A to B. The flux of link types is then reversed, and 
the steady state network structure is similar (but not 
identical) to that expected from the community structure 
parameters d, /. The infection persists at high levels in 
both communities at steady state. 

We can estimate the time until infection incursion in the 
B community as follows. From Equations (|3I6I10I) . the 

fraction of AA links P aa = PsaSa + Psala + Piaia evolves 
according to 



1 [aPsalb ~ (1 - a)PsaIa] 



(13) 



Since we are interested in the critical time when infection 
starts to spread in community B, we can assume Pjb ~ 
and hence Psaib ~ 0. Thus 



-W(l - a)P Sa Ia 

-w(l - a)<yP aa 



(14) 



where 7 = Psaia/Paa- Since / is close to zero, we can 
use the single mean field equations to approximate 7, the 
fraction of SI links in community A (see Appendix A). 
We can then solve Equation (fT4"|) for the critical time t c 
for infection incursion if we know P aa at that time. Since 
we can predict the critical average degree for community 
B in order for the disease to spread, we can also find P aa 
at that point. However, for the full system, the infection 
in community B starts to spread much earlier than the 
time found by using mean field equations because of the 
stochastic nature of our model. Even so, we can solve 
the equation for t as follows: 



t = In 



Paa{0) 



Paa(t) 



[w(l - a)j}~ 



(15) 
(16) 



from which it can be shown that the critical time t c is 
proportional to 4 . In Figure [3 we can see that the rela- 
tionship t c oc j holds for the full system. 
When the system reaches steady state (late time in Fig- 
ure U) , the network geometry does not return exactly to 
that expected from the community structure parameters 
/, d because the SI links are not uniformly distributed. 
This effect is most pronounced when the infection rate is 
below the critical infection rate for the low connectivity 
network, because then the SI link distribution is the most 
nonuniform. Figure |5] shows the steady state commu- 
nity structure versus the infection rate. Deviations from 
the community structure specified by /, d increase as the 



6 



**- 10 





x 10 



FIG. 5: The critical time for infection incursion in community 
B vs fraction of cross links /. The critical time values are 
averaged over 5 runs except for / = 1CP 6 . In each run, the 
time when the fraction of infecteds in community B exceeded 
0.1% is considered the critical time. 



FIG. 6: Fraction of links vs infection rate p. Curves: mean 
field steady states; points: Monte Carlo steady state averages. 
Black: fraction of AA type links; gray: fraction of BB type 
links. (The fraction of AB links is less than 1%.) N a = 5x 10 3 , 
N b = 5 x 10 3 , K = 10 6 , / = 1(T 4 , d = 0.7, w = 0.04, 
r = 0.002. The data is from the simulations done in Figure [3] 



infection rate approaches the Hopf bifurcation point at 
p w 0.0019. Thus the steady state average degree ob- 
served in the two communities in the presence of an epi- 
demic can be different than that expected in the absence 
of an epidemic. The adaptation has a homogenizing ef- 
fect, bringing the degrees in the communities closer to 
each other. 

For d = 0.9, the behavior of the system is very interest- 
ing. In Fig. [3 we started with similar initial conditions 
as for d = 0.7 in Fig.HJ The system behaved similarly for 
a long time, but once the infection level reached a high 
value in community B, it could not stay in that state, 
because a higher d means a higher out flux rate from BB 
type links. The average degree in B started to decrease, 
causing infection to die out again in B. A periodic so- 
lution with a long period is observed for a range of p 
values. 



IV. CONCLUSIONS 

We have studied epidemic spread in a network of two 
communities with different average degrees. Cases with 
and without disease avoidance rewiring were considered. 
Rewiring rules were chosen so that the community struc- 
ture would be preserved if links between susccptibles and 
infectives occurred uniformly throughout the network. 
The steady state bifurcation structure was obtained for 
static and adaptive cases. A mean field theory based on a 
moment closure approximation accurately predicted the 
steady state infection levels and network structure ob- 
served in stochastic simulations of the full model. 
In the static network case, weakly connected communi- 
ties displayed significantly different infection levels. Low 
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FIG. 7: Time series from Monte Carlo simulation of full sys- 
tem, (a) Infected node fractions. Black curve: fraction of 
infecteds in community A; gray curve: fraction of infecteds 
in community B. (b) Fraction of links. Black curve: frac- 
tion of AA type links; gray curve: fraction of BB type links. 
Initial condition is as described in the text. N a = 5 x 10 3 , 
N b = 5 x 10 3 , K = 10 6 , / = 10" 4 , d = 0.9, w = 0.04, 
r = 0.002 and p = 0.0014. 



infection levels could persist in a subthreshold commu- 
nity weakly connected to a high degree, high infection 
community. Increasing the number of connections be- 
tween communities led to more similar behavior of the 
two communities. In contrast, communities in adaptive 
networks displayed similar infection levels even if weakly 
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connected. Steady states with high infection in one com- 
munity and low in the other did not exist for the adaptive 
network case. 

The absence of steady states with significantly differ- 
ent infection levels was explained by considering network 
adaptation in the presence of nonuniformly distributed 
SI links. If one community has few infectives, there is a 
net flux of links into that community until its degree is 
high enough to support the infection. We estimated the 
time until this infection incursion based on mean field ar- 
guments and found that the time increases as the number 
of cross links between communities decreases. 
We also observed changes in the steady state network 
geometry due to adaptation in the presence of infec- 
tion. These changes were most significant near bifurca- 
tion points. The adaptation tended to bring the average 
degrees of the communities closer to each other. Thus the 
adaptation promotes greater similarity between commu- 
nities in both network structure and infection levels. 
The model presented in this paper is the first to include 
community structure in epidemic spread on an adaptive 
network. Future work is needed to extend the model 
to more realistic scenarios. For example, the number of 
communities could be increased beyond two. We have ob- 
served that the convergence time to steady state can be 
very long for weakly coupled communities, so it is possi- 
ble that an epidemic would not reach steady state during 
physically realistic time scales. Thus the transient be- 
havior should also be studied in more detail. Identifying 
when communities become at risk for incursion of infec- 
tion could be valuable in knowing when epidemic control 
measures are needed. Another area for future extension 
is to change the rules for adaptation, such as cutting or 
temporarily deactivating links rather than rewiring them. 
This work was supported by the Army Research Office, 
Air Force Office of Scientific Research, and by Award 
Number R01GM090204 from the National Institute Of 
General Medical Sciences. The content is solely the re- 
sponsibility of the authors and does not necessarily repre- 
sent the official views of the National Institute of General 
Medical Sciences or the National Institutes of Health. 



Appendix A: Analytical solution for a single 
community 

For a single community (/ — 0), the mean field equations 
are: 



Pi 



pK 

-TP; + P W P SI 



P SI = 2r{l-P SI -P ss )-{r+p + w)P S i 

2pK PssPsi pK PIj 

N P s N P s 

Pss = rP S i + wP S i 
2 P K P ss Psi 

N P s 
At steady state, we obtain 

-r Pl + fPs, = 



2r(l - Psi - Pss) -(r+ P + w)P : 



si 



2pK PssPsi pK P§j _ Q 



rP : 



N P s N P s 

2pK PssPsi 



si 



wPsi - 



N Ps 



= 



(Al) 

(A2) 
(A3) 



To find the endemic steady state, we first solve |A1| for 
Psi, and substitute into |A2| . Then we solve |A2| for 
Pss in terms of Pi. After substituting that into |A3| . wc 
obtain a quadratic equation in Pi, APf + BPi + C = 0, 
where 



A 
B 

C 



p — w 

2w -p- 2p^ 



K 

N 



w 



The quadratic can be solved analytically for Pi, and then 
Psi and the other link variables can be computed. 
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